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1. Introduction 

The method of cluster expansions was introduced in the 1930's in statistical me- 
chanics in order to study gases of classical interacting particles. Its main achieve- 
ment, from the point of view of physics, may be the derivation of the van der Waals 
equation of state for weakly interacting systems. The method was made rigorous 
by mathematical-physicists in the 1960's, see [23] and references therein. 

The method split afterwards. One branch involves continuous systems, with ap- 
plications to classical systems [THl HZl [S] , quantum systems [TUl El [21] , or quantum 
field theory [T^l [IHl [H [S] • The other branch involves polymer systems, i.e. discrete 
systems with additional internal structure [131 [T] [31 [TSl [13 [S] [S] ■ An important 
step forward was the article of Kotecky and Preiss with its simplified setting and 
its elegant condition for the convergence of the cluster expansion [15j . 

The methods for proving the convergence are diverse. Let us mention the study 
of Kirkwood-Salsburg equations that involves correlation functions, see [23] and 
references therein; the algebraic approach of Ruelle [23 ; combinatorial approaches 
using tree identities [THl [3 [21 [S] ; inductions for discrete systems [7| |31 HH] ■ 

Important and useful surveys were written by Brydges [4] , Pfister [20] , Abdesse- 
1am and Rivasseau [T] . Recent articles have been devoted to combinatorial aspects 
[251 [HI m] and to weakening the assumptions [S] [5] [22] ■ 

The method of cluster expansions applies when the objects do not interact much; 
this the case when they are far apart (low density), or when interactions are weak. 
An extension of the criterion of [15] that takes into account these two aspects was 
proposed in [26]; it apphes to both discrete and continuous systems. 

All abstract (i.e. general) approaches involve restrictions that correspond to re- 
pulsive interactions. Yet the old results for classical and quantum systems only 
assume stable interactions, that may include an attractive part. The aim of the 
present article is to propose a general approach that applies to discrete and contin- 
uous systems with repulsive or stable interactions. Our proof is split into several 
independent steps and this helps clarify the situation. 

The setting and the results are presented in Section [21 We consider applications 
to classical systems of particles in Section [31 to polymer systems in Section [H and 
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to the quantum gas in Section [5l A fundamental tree estimate is derived in Section 
[6l and the theorems of Section [2] are proved in Section [71 

2. Cluster expansions 

We consider a set X whose elements may represent widely different objects — in 
the three applications considered in this article, an element a; £ X represents (i) the 
position of a classical particle, (ii) a polymer, i.e. a connected set of Z'', and (iii) a 
closed Brownian bridge. For the general abstract theory, we assume the structure 
of a measure space, {lL,X,fi), with fi a complex measure. We denote \fi\ the total 
variation (absolute value) of ^. Let u and C be complex measurable symmetric 
functions on X x X, that are related by the equation 

C(x,2;)= e-"(-'^) -1. (2.1) 

We allow the real part of u to take the value +oo, in which case Ci^iU) = ^1- In 
typical applications u{x,y) represents the interactions between x and y, and the 
value +00 corresponds to a hard-core repulsion. We define the "partition function" 

by 

Z = ^ — dfi{xi)... / d^(a;„)exp|- ^ M(xi,a;j)|, (2.2) 

n > 1 ^ i<j ^ Ti 

or, equivalently, 

M^i)--- j d^a^n) n (1 + C(a;.,xj)). (2.3) 

The term n = of the sums is understood to be 1. 

The main goal of cluster expansions is to express the partition function as the 
exponential of a convergent series of "cluster terms" . The main difficulty is to prove 
the convergence. We first assume that the potential u is stable. 

Assumption 1. There exists a nonnegative function b on \ such that, for all n 
and almost all xi, . . . ,Xn G X, 

n 

1 ^ i<j ^ n i—1 

In other words, we assume the lower bound 

n 

Reu{x,,Xj) ^ ~Yb{xi). (2.4) 

1 ^ i<j ^ n i—1 

When the function b is constant, this is the usual definition of stability. "Almost 
all" means that, for given n, the set of points where the condition fails has measure 
zero with respect to the product measure If X is countable, the condition 

must be satisfied for all cci, . . . , a;„ such that iJ,{xi) ^ 0. 

The second condition deals with the strength of interactions. 

Assumption 2. There exists a nonnegative function a on H such that for almost 
all X e X, 

d\^l\{y)\ax,y)\e^<^y^+''^(y'> < a{x). 
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In order to guess the correct form of a, one should consider the left side of the 
equation above with a{y) = 0. The integral may depend on x; a typical situation 
is that X is characterized by a length i{x), which is a positive number, so that the 
left side is roughly proportional to £{x). This suggests to try a{x) = c£{x), and one 
can then optimize on the value of c. 

We also consider an alternate criterion that involves u rather than C. It is inspired 
by the recent work of Procacci [22] . Let 

-/ ^ juix^y) i{Reuix,y) ^ CO, 

u(x,y) = < . (2.5) 

I 1 if Re u{x, y) — oo. 

Assumption 2'. There exists a nonnegative function a on X such that for almost 
all X £lL, 

d|/i|(y)Kx,2/)|e'^(^)+^(^) s; a{x). 

For positive u we can take b{x) = 0; and since 1 — e^" ^ u, Assumption [2] 
is always better than Assumption [2l We actually conjecture that, together with 
Assumption [1] a sufficient condition is 

d\ti\iy)min{\ax,y)l\u{x,y)\)e^^y^+'^y^ sCa(x). (2.6) 

That is, it should be possible to combine the best of both assumptions. In this 
respect Assumption [2] is optimal in the case of positive potentials, and Assumption 
[2lis optimal in the case of hard core plus negative potentials. 

We denote by Qn the set of all graphs with n vertices (unoriented, no loops) and 
Cn C Gn the set of connected graphs with n vertices. We introduce the following 
combinatorial function on finite sequences {xi, . . . , x„) of elements of X: 

I 1 if n = 1, 

^(xi,...,x„) = <^ (2.7) 

The product is over edges of G. 

Theorem 2.1 (Cluster expansions). Suppose that Assumptions]^ and\^ or{l\and 

[2\ hold true. We also suppose that J d|/^|(?/)| e"^^^"'"^''^'') < oo. Then we have 

Z = exp| ^ d^j.{xi)...dfi{xn)(p{xi,...,Xn)y (2.8) 

n ^ 1 

The term in the exponential converges absolutely. Furthermore, for almost all Xi G 
X, we have the following estimate 

J2nJ d\fi\ix2) . . . J d\^,\ixn)Hxu■■■,Xn)\ < ( c^^^^^ - 1) e^'^^^K (2.9) 

n ^ 2 

(Under Assumption [13 Eq. holds with e''^^^) instead of e^^^^i) .) 

Let us turn to correlation functions. We only consider one-point and two-point 
correlation functions since these are the most useful and expressions become more 
transparent. We refer to [IS] for more general functions. First, we define the 
unnormalized one-point correlation function by 

Z{xi)=Y.^^^^ljM^^)---jAKxn) n (l + C(2;.,a^,)) (2.10) 



1 ^ i<3 ^ " 
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(the term n = 1 is 1 by definition). And we define the unnormalized two-point 
correlation function by 

Z(X1,X2)= in-2-)\ I dM(^3)--- / n {l + Q{^r,Xj)) (2.11) 

(the term n = 2 is equal to 1 + ({xi,X2))- Notice that Z{xi) can be viewed 
as a regular partition function, given by Eq. (|2.2p . but with the modified measure 
(l+^(a;i, x))^{x) instead oi ^{x). The normalized correlation functions are Z{xi)/Z 
and Z{xi,X2) / Z . As is shown in Theorem 12.21 they can be expressed using the 
"cluster functions" 

Z{xi) = X! " / '^l^^^'i) ■ ■ ■ j ^t^i^n) f{xi, . . . , Xn), (2-12) 

and 

Z{xi,X2)^ ^n(n-l) dfi{x3) . . . dfi{xn)^p{xi, . . . ,Xn)- (2-13) 



n>2 



Notice that I Z(xi) I Qa(x^)+2b(x^) by ([M 



Theorem 2.2 (Correlation functions). Under the same assumptions as in Theorem 
\2.1{ we have 

-- Z(xi), 



Z 

Z{xi,X2) 

z 



= Z{xi)Z{x2) + Z{xi,X2)- 



In statistical mechanics, the relevant expression is the truncated two-point cor- 
relation function 

Z{xi,X2) Z{xi)Z{x2) 
Z Z2 ■ 

When the cluster expansion converges, it is equal to Z[xi , ^2) by the theorem above. 
This function usually provides an order parameter for phase transitions and it is 
useful to estimate its decay properties. 

Theorem 2.3 (Decay of correlations). If AssumptionsUi and[B hold true, we have 
for almost all x,y 



\Z{x,y)\ < e°(^)+2''(2') [|C(x,?/)|e°(^)+2«'(^ 

d\fi\ix,)... dH(x„)niC(2:.,x,+i)|e'^(- 



m 

)+2b(xi) 



m> 1' 



(with xq = X and Xm+i =y)- If Assumptions[I\ and\2} hold true, we have the same 
bound but with \u(-, ■)] instead of ■)], and e''^'^ instead of e^^''^ . 

In many applications the functions Cix, y) and u{x, y) depend on the difference 
X — y (this assumes that X has additional structure, namely that of a group). The 
estimates for \Z(x,y)\ are given by convolutions. 

The theorems of this section are proved in Section [71 
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3. Classical gas 



We consider a gas of point particles that interact with a pair potential. We 
work in the grand-canonical ensemble where the parameters are the fugacity z and 
the inverse temperatm'e /3 (both are real and positive numbers). The set X is an 
open bounded subset of R*^ and [i{x) = zdx with da; the Lebesgue measure. We 
actually write A = X so as to have more traditional notation. The interaction is 
given by a function f7 : R*^ ^ M U {00} which we take to be piecewise continuous; 
u{x,y) = f3U{x — y). We suppose that U is stable, i.e. that there exists a constant 
B ^ such that for any n and any xi, . 

f/(x, 



E 

1 ^ i<j ^ n 



Xj) ^ —Bn. 



(3.1) 



Our Assumption [T] holds with b{x) = (3B. The system is translation invariant so all 
cc e M'' are equivalent. The function of Assumptions [2l and W\ can then be taken to 
be a constant, a{x) = a. We seek a condition that does not depend on the size of 
the system. Then integrals over y are on R'' instead of A. By translation invariance 
we can take x — 0. 

Assumption [2] gives the condition 



ze 



2f3B 



-0U(v) 



1 dy ^ a e 



(3.2) 



We obviously choose the constant a that maximizes the right side, which is a = 1. 
This condition is the one in :23] . Let us now assume that U consists of a hard core 
of radius r and that it is otherwise integrable. Again with a = 1, Assumption [2l 
gives the condition 



ze 



s; e" 



(3.3) 



P \Uiy)\dy 

J\y\>r J 

Here, |B| = 7r^/'^/T{^ + 1) is the volume of the ball in d dimensions. This condition 
is often better than (|3.2p . Without hard core it is the one in [S]. The domains 
of parameters where these conditions hold correspond to low fugacities and high 
temperatures. 

The thermodynamic pressure is defined as the infinite volume limit of 

PAiP,z) = j^logZ. (3.4) 



|A| 



Using Theorem I2.1[ we have 
1 



|A| 



dxi 



■n > 1 



dxy 



dXn(p{xi, ...,Xn) 



(3.5) 



Consider now any sequence of increasing domains Ai C A2 C . . . such that A„ 
R''. Thanks to the estimate (|2.9[) . and using translation invariance, we get 



p{(3, z) = lim pA„ (/3, z) = V z" / dx2... 



dXnipiO,X2,..-,Xn)- (3.6) 



(The term with n = 1 is equal to z.) This expression for the infinite volume pressure 
p{P, z) should be viewed as a convergent series of analytic functions of /3, z. Then 
p{(3, z) is analytic in (3 and z by Vitali theorem and no phase transition takes place 
in the domain of parameters where the cluster expansion is convergent. 
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The truncated two-point correlation function (t{x) is given by Z{0,x). We con- 
sider only the case of Assumption [2] but a similar claim can be obtained with 
Assumption [21 Let c{x) be a function that satisfies the triangle inequality. The 
estimate of Theorem 12.31 yields 

e<^K{x) e^+^/ss ^'"e'"+2'^^"'(e^(-)|e-''^(-) (3.7) 

m ^ 

(with r° = /). Recall that ^ ||/||oo||/|ir \ and let 

||e=(-)|e-'3^(-) -1|||^. (3.8) 

Then we get 

aix) ^ e-^(^) e2+4^^S^(l-^ei+2'3^Ci)"\ (3.9) 

If the inequality (|3.2p is strict, one can usually find a function c{x) that satisfies 
the triangle inequahty and such that Ci < {ze^~^^^^ )~^; the truncated two-point 
correlation function then decays faster than e"^^^-* . 



4. Polymer systems 

Polymer systems are discrete, which is technically simpler, but they also have 
internal structure. The first application of cluster expansions to polymer systems 
is due to Gruber and Kunz [13] . Among the many articles devoted to this subject, 
let us mention [THl El [S] ■ The main goal of this section is to illustrate our setting; 
we therefore restrict ourselves to a specific model of polymers with both repulsive 
and attractive interactions. 

Our space X is the set of all finite connected subsets of Z'*. The measure fj, is the 
counting measure multiplied by the activity z{x) (a function X ^+)- We choose 
z{x) = e~'''l^l with 7 > 0. The interaction is hard core when polymers overlap and 
it is attractive when they touch: 

joo if xny ^ 0, 

u{x,y)^i ^ ^ -f n (h (^-^^ 
\—ric(x,y) it x n y = It). 

Here, c{x, y) is the number of "contacts" between x and y, i.e. the number of bonds 
between sites of x and y; > is a parameter. The interaction is zero when the 
distance between polymers is greater than 1. 
The stability condition can be written 

n n 
i—l j^i i—1 

Only disjoint polymers need to be considered, the left side is infinite otherwise. The 
sum over j is always larger than —77 times the number of bonds connecting Xj with 
its exterior. Thus we can take b{x) — r]d\x\. 

The function a in Assumption [2l grows like |a;|, so it is natural to choose a{x) = 
a\x\ for some constant a. A sufficient condition is that 

z(y)e'^l«l+'"*l^l + J2 ?7z(y)c(a;,?/)e°l^l+'"^l''l a\x\. (4.3) 

y^yr\x^^ y,dist {x,y) — l 
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We can bound ric{x,y) by 2rjd\y\. Summing over the sites of x, and requiring that 
y contains the given site or comes at distance 1, we get 

^(l + 2(i??|y|)z(2/)e°l''l+'"*l^l a. (4.4) 

We used the fact that the activity is translation invariant. If a; is a connected set, 
there exists a closed walk with nearest neighbor jumps whose support is a;, and 
whose length is at most 2(|x| — 1). This can be seen by induction: knowing the 
walk for a;, it is easy to construct one for x U {«}. The number of connected sets of 
cardinality n that contain the origin is therefore smaller than the number of walks 
of length 2n — 3 starting at the origin, which is equal to (2(i)^"~^. Then it suffices 
that 

g-„(7-a-3d„-21og2d) ^ (2dfa. (4.5) 

n ^ 1 

This is equivalent to 

7 ^ a + log(l + (2^) + 3d?7 + 2 log 2d. (4.6) 

Assumption W\ holds for any a. Using log(l ~\-t)^t and optimizing on a, we find 
the sufficient condition 

7 > 2 (2d) ^3/2 + 3d?? + 2 log 2d (4.7) 

with a = (2d)-3/2. 

We have just established the existence of a low density phase provided the activ- 
ity is small enough. The condition depends on the contact parameter rj. For large 
rj one should expect interesting phases with many contacts between the polymers. 



5. Quantum gas 

We follow a course that is similar to Ginibre [11] , using the Feynman-Kac formula 
so as to get a gas of winding Brownian loops. We get comparable results, with a 
larger domain of convergence in the case of integrable potentials. Winding Brownian 
loops are kind of continuous polymers; they combine the difficulties of both cases 
above — the continuous nature and the internal structure. 



5.1. Feynman-Kac representation. The state space for TV fermions (resp. bosons) 
in a domain A C M'' is the Hilbert space L^„ti(A^) (resp. Lly^{K^)) of square- 
integrable complex functions that are antisymmetric (resp. symmetric) with respect 
to their arguments. The Hamiltonian is 

N 

H = -J2^.+ E U{q,-q,), (5.1) 

i=l 1 ^ i<j ^ N 

with Ai the Laplacian for the i-th variable and U{q) a multiplication operator. As 
in the classical case, we consider the grand-canonical ensemble whose parameters 
are the fugacity z and the inverse temperature (3. The partition function is 

Z= J2 ^"^Tr e-f^" . (5.2) 
JV ^ 
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We need to cast the partition function in the form (12. 2p , which can be done using 
the Feynman-Kac representation. Namely, we have 



N ^0 ■ TrG5jv J A J J 

(n^A(a;.))exp{-i J2 / C/(w,;(s)-w,(s))ds}. (5.3) 



iSat is the symmetric group of N elements; £{n) is equal to the signature of the 
permutation tt for fermions, £(7r) = 1 for bosons; VF* ^/ is the Wiener measure for 
the Brownian bridge from q to q' in time t — the normalization is chosen so that 

dWl^,{Lu) = (27rt)-'^/2c^l9-«'l'/2* ; (5.4) 

Xa('^) is one if w(s) G A for all ^ s ^ 2/3, it is zero otherwise. An introduction 
to the Feynman-Kac formula in this context can be found in the survey of Ginibre 

m- 

The right side of Eq. (|5.3p is well defined for a large class of functions U :M.'^ 
M U {oo}, that includes all piecewise continuous functions. Thus we take (|5.3p as 
the definition for Z. Under additional assumptions on U, (|5.3p is equal to (|5.2p 
with Hamiltonian (|5.ip and with Dirichlet boundary conditions. 

We now rewrite the grand-canonical partition function in terms of winding loops. 
Let Xfc be the set of continuous paths [0, 2pk] — > M*^ that are closed. Its elements 
are denoted x — {q,k,uj), with q G M'' the starting point, k the winding number, 
and Lo the path; we have lo{0) = Lj{2j3k) = q. We consider the measure /i given by 

/c+l k 

^{^Ax) = L_^d<ze-"(-' XA(c.Xf (d^). (5.5) 



Here, v{x) is a self-interaction term that is defined below in (|5.7p : e = — 1 for 
fermions and -f 1 for bosons. Let X = ^ iX^; the measure ^ naturally extends 
to a measure on X. The grand-canonical partition function can then be written as 

n 5s ' 1 ^ i<j ^ n 

Let X — {q,k,uj) and x' — {q',k',uj'). The self-interaction v{x) and the 2-loop 
interaction u{x, x') are given by 

.2/3 



r,(x) = ^ i/ [/(t^(s + 2/3^) -cj(s + 2/3m))ds 

fc-lfc'-l „2/3 

u(a;, a;') = E E 5 / ^("^(^ + ^Z?^) " + W))ds. 



We now treat separately the case of integrable potentials and the case of more 
general potentials. 
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5.2. Stable integrable potentials. We suppose that U is stable witli constant 
B, i.e. it satisfies Eq. (13. For given loops xi, . . . ,a:„, stability implies that 



^v{xi)+ u{x,,Xj) ^ -pB^k,. (5.8) 

2—1 1 i<j ^ n 2—1 

Then Assumption [1] holds with 

h(x) ^ (3Bk + v{x). (5.9) 

Notice that b{x) ^ 0, again by the stability of U . 

We use Assumption W\ and we choose a function a{x) — ak with a constant a to 
be determined later. Explicitly, the assumption is that for any x = {q, k,Lu) 



■^(■'■'')+f>('-^")\u{x,x')\ ak. (5.10) 



We have lifted the restriction that uj'{s) G A because we want a condition that does 
not depend on A. Eq. (|5.10p is easier to handle than its appearance suggests. Notice 
that the term in the second exponential is just j3Bk' . Using J dH^^^(u;)/(w) = 
/ dWoo(t^)/(t^ + ?), and the definition (j5.7|) of u{x,x'), it is enough that 

1=0 -^^ k' ^ 1 r=o * 

^2k'l3 



dWo^a^'{Lu')\U{uj{s + 2p£)-uj\s + 2p£')-q')\ ak. (5.11) 

We can immediately integrate over q' , which yields ||C/||i. The Wiener integral then 
gives (47rfc'/3)~''/^ and we get the equivalent condition 



k' {a+l3B)k' 



For any d, /?, ||C/||i, a, the inequality holds for z small enough. Notice that z < 1 in 
any case. One can get a more explicit condition for d ^ 3 by choosing a such that 

niiC(f) 



r flfl^ilnMlJ ^ nil rt^1•?^ 



(477/3)''/ 

Here, ({^) — J^n > ^^e Riemann zeta function. 

When d — 3 and when the potential is repulsive, one can rewrite (|5.13p in a 
more transparent way. Let — ^||f7||i denote the Born approximation to the 
scattering length. The condition is then 

z ^ exp{-^^}. (5.14) 

The critical fugacity is expected to be greater than 1. The present result helps 
nonetheless to obtain a range of densities where the pressure is analytic. In the 
bosonic case it compares well with physicists' expectations 
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5.3. Stable potentials with hard core. The presence of a hard core makes the 
situation more complicated; we only sketch the argument in this section without 
trying to get explicit bounds. Our aim is to show that, using Theorem 12. 1[ the 
problem of convergence of the cluster expansion reduces to estimates of Wiener 
sausages. 

We consider an interaction U = U' + U" . We assume that U' > is a repulsive 
potential of radius r, with a hard core of radius < r' ^ r, that U"{q) = 
for \q\ < r, and that U" is integrable otherwise. We suppose that the stability 
condition takes the slightly stronger form 

n 

J2u{qo-qi) ^ -B (5.15) 

1=1 

for any go, . . . , g„ such that \qi — qj\ > r' . For potentials with a hard core this is 
equivalent to the property (|3.ip . possibly with a different constant B. Then one 
has [lO] 

n n 

^v{x^) + 2 ^ u{xi,x^) ^ -2(3B^h. (5.16) 

i=l 1 ^ i<j ^ n j=l 

Then Assumption [T] holds with 

b{x) ^ PBk+ ^v{x). (5.17) 

Notice that the stability condition also holds with b given in (|5.9|) (and with a better 
constant B). The advantage of (|5.17p is the factor i in front of v{x). Then e^*"^^^ 
involves a term that cancels the self-interactions of x. 

Given a loop x — {q,k,uj), let S{x) be the Wiener sausage generated by a ball 
of radius r when its center moves along the trajectory uj: 

Six) = {qeW^ : \uj{s)-q\ < r for some s £ [0,2fc/3]}. (5.18) 

We also define the Wiener sausage that corresponds to the difference of two winding 
loops X — {q, fc, uj) and x' — (q' , fc', uj'): 

S{x,x')^ {q" ■.\uj{s + 2f3£)~uj'{s + 2l3£')-q"\ s$ r 

for some £ = 0, ... ,k -!;£' = 0, ... ,k' - I; s € [0,2(3]}. (5.19) 
We denote the volume of a Wiener sausage S{-) by \S{-)\. One can check that 

with |B| the volume of the unit ball (see Appendix 2 in [lO]). 

We choose a{x) = \S{x)\ in Assumption[21 Then a sufficient condition is that 
for any a; G X, 

^ / d,'/d<r(-')IC(x,-' + 9')|el"(^)l ^\S{x)\+k. 

(5.21) 

We consider separately the cases where q' belongs or not to S{x,x'). First, 

\ax,x' + q')\dq' ^ \S{x,x')\, (5.22) 
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which we bound using ()5.20p . Second, using | e ) — 1| ^ e'^^'^ , 

\Cix,x' + q')\dq' ^ l3\\U"\\ikk'e^^''' . (5.23) 



lR<'\S{x,x') 

We certainly get (|5.2ip if we have the two incquahties 

yk' (2/3B+1)*:' 



E /d</(c.')el-(^')l|5(x')| =^ 

:' 5s 1 

E ^''^'e(^^^+^)'='/d<f(c.')el^(^')l ^ 



(5.24) 



fe' > 1 



One can estimate the integrals of Wiener sausages, see [TU], so that both conditions 
hold if z is small enough. 

Now that the cluster expansion is known to converge, it is possible to write the 
pressure as an absolutely convergent series of analytic functions in /3 and z. It is 
also possible to study the decay of correlation functions. In the case of potentials 
that consist of hard core plus integrable part, one can apply Assumption 2' instead. 
This may give better results, especially if the integrable part is mostly attractive. 

6. Tree estimates 

In this section we obtain estimates of sums of connected graphs in terms of sums 
of trees. Our main result is Proposition 16.11 below. Such estimates seem to have 
been introduced by Penrose [19] and they have often been considered in the past 
[51 [21 m ini dOl HI IH] • We introduce a minimal setting that clarifies its role in the 
cluster expansion. Namely, we fix the polymers so we only deal with the numbers 
that represent their interactions, C or m, and the stability function b. Assumption 
[T]is vital here, but Assumptions and [2l are not used in this section. 

Let Tn G Cn denote the set of trees with n vertices. Let n be an integer, &i, . . . , 6„ 
be real nonnegative numbers, and Qj = Qi, I ^ i,j ^ n, be complex numbers. 
We assume that the following bound holds for any subset / C {1, . . . , n}: 

n ^ He'-: (6-1) 

,i<j i€l 

Let Uij be such that Qj — e^"'^ — 1; let Uij = 1 if = — 1, and Uij — Uij 
otherwise. We state two distinct tree estimates, the first one involving \ and the 
second one involving \uij \ . These bounds will allow to prove the convergence under 
either Assumption [2] or Assumption \T\ 

Proposition 6.1. // (|6.ip holds true, we have the two bounds 

n 

(-) |E n c.| ^ (n^"*) E n ici; 

n 

(b) |E n c.| ^ (n^^o E n i^^.i- 

Gec„ {jj}eG 1=1 GeT„ {i,j}eG 

We actually conjecture that the following estimate holds under the same hy- 
potheses: 

n 

|E n ^ (n^'*) E n niin(Ki,ic.,i). (6.2) 

Gec„ {i,i}eG i=i GGr„ {ij}eG 
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We prove Proposition 16. II (a) below using Ruelle's algebraic approach 23\. This 
method is usually combined with a Banach fixed point argument for correlation 
functions. However, we use it differently so as to get a tree estimate. Proposition 
16.11 (h) follows from a tree identity due to Brydges, Battle, and Federbush [5l[2j|4], 
combined with an argument due to Procacci |22| ; its proof can be found at the end 
of this section. 

Let A be the set of complex functions on the power set V{{1, . . . ,n}). We 
introduce the following multiplication operation for f,gGA: 

f^giI)^Y.f^^^9{I\J)- (6.3) 

JC/ 

We use the standard conventions for sums and products, namely that the empty 
sum is zero and the empty product is one. Together with the addition, ^ is a 
commutative algebra with unit 1a{I) ~ (5/ g. It is possible to check that each 
f £ A \ {0} has a unique inverse, which we denote We have 

r\i)- E f{Ji)---fiJk). (6.4) 

Jl,...JkCl 

JinJj=0,uJi=7 

Let Ao be the subset of functions / such that /(0) — (^o is an ideal of A). Notice 
that f *'^ = for any k > n , when f £ Aq. We define the exponential mapping 
exp_4 : Ao ^ Ao + 1a by 

exp^/ = l^ + /+ir' + --- + ^r". (6.5) 

Let $ and ^ be the functions defined by 

E n c.,' 

GeC(7) {Lj}eG 

n (i+c.,)= E n c.,. 

i,J&I,i<J G<£Q{I) {i,j}i£G 

Here, G{I) (resp. C(/)) is the set of graphs (resp. connected graphs) on /. We have 
the relation 

* = exp_4 $. (6.7) 
We also introduce an operation that is reminiscent of differentiation: 

D.,m J (6,S) 
I otherwise. 

One can check that D^^y cxp^ f ~ (exp_4 /) * 
For disjoint /, J C {1, . . . , n}, we define 

5(/, J) *i?/*)( J). (6.9) 

Let J C {1, . . . , n}. The assumption of Proposition 16. II implies that 

n n ^ n^""- (^.lo) 

iei jei\{i} iei 
Then there exists i G I such that 

n + ^ e^"-. (6.11) 
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Such i is not unique in general but it does not matter. We consider a function l that 
assigns one of the indices i — above to each nonempty subset / C {1, ■ • ■ ,n}. 
Notice that G / for any subset /. It is also useful to introduce the notation 
/' = /\ {.(/)}. 

Lemma 6.2. The function g {I , J) of Eq. (|6.9p is solution of the following equation. 

7^,^)= (n(i+W))) E(n^M(/))5(/'uif,j\i^) *//^0. 

iei' KcJ ieK 

Since the equation gives J) in terms of g{K, L) with \K\ + \L\ = \I\ + \ J\ — 1, 
it is well defined inductively and it has a unique solution. Notice that .g(0, 0) = 1, 
and that g{{i}, 0) = 1 for any index i. 

Proof. Recall the definition (|6.6p of 5*. For disjoint /, K we have 

vi/(/uic)=( n (1+ W)))*(^'u^) 

(6 12) 

= (n(i + W)))(E \{CK.ii))w^K). 

jel' LCK keL 

Then 

g{I,J) = ^'*'--'^\j\K)'^{lUK) 

KCJ 

= (n(l + W))) E {l['^^Mi))'^*^''\j\K)^iI'^K) 

jel' LCA'CJ keL 

n(l + W))) Edl^M/)) E ^*^~'HJ\L\K')^{I'ULUK'). 

jel' LCJ keL K'GJ\L 

(6.13) 

The last sum is equal to g{I' U L, J \ L). One recognizes the equation of Lemma 
iH □ 

We now estimate the function g using another function h that satisfies an equa- 
tion that is similar to that of Lemma 16.21 



'hi9,J) = 5^j, 

h{I,J)= e^'^in J2(^Y[\Qai)\)HI'^K,J\K) if/^0. (6-14) 

KCJ iCK 

It also has a unique solution. Since Hie/' |1 + Ci.t(/)l ^ e^**'!^) , we can check 
inductively that 

\g{I,J)\ «C h{I,J) (6.15) 
for any sets /, J (with I D J — 0). Now the function h can be written explicitly 
[T71 W] . Let J^i[J) be the set of forests on / U J rooted in I. That is, a graph 
G S ^i{J) is a forest such that each tree contains exactly one element of /. 

Lemma 6.3. The solution of Eq. (|6.14p is 

M/,j)=( n E n icd- 

icivjJ GeJ^iiJ) {i,j}eG 
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Proof. Since the solution to Eq. ()6.14p is unique, it is enough to check that the 
Ansatz of the lemma satisfies the equation. First, let us observe that both sides are 
multiplied by HiG/u/ ■ Thus it is enough to consider the case bi = 0. 

The sum over graphs in Ti{J) can be realized by first summing over the set K 
of indices (necessarily in J) that are connected to i(/); then over sets of trees in 
J\ K, and over connections to /' U K. Explicitly, 

E nicd- E(niW)i) E n ic^- (6.16) 

G£ri{J}{i,j} KCJ i<£K G<£ri,uKiJ\K){ij}€G 

This equation is precisely (|6.14p . □ 

Proof of Proposition \6.1\ (a). When / has a single element, the function g is equal 
to 

5({1},{2,...,7^}) Cy- (6.17) 

Gec„{i,j}eG 

This is the left side of Proposition 16.11 (a). We have J-{i^{{2^ . . . ,n}) — Tn, the 
set of trees with n vertices. Thus ^({1}, {2, . . . ,n}) is equal to the right side of 
Proposition [nH] (a), and the proof follows from Eq. (|6.15p . □ 

We now turn to the proof of Proposition l6.1l (b). Notice that in absence of "hard 
cores", i.e. when dj ^ —1, our claim is just a reformulation of Corollary 3.2 (a) of 
[4|. The present proof follows [22] . 

Proof of Provosition \6.1\ (h). Let P be the set of {i, j} such that Qj — —1, i.e. such 
that Reuij — oo. We regularize those numbers by setting 

\Uij lf{l,j}^P. 

This allows to use the tree identity of [51 H] ; we will eventually take m to infinity. 
The tree identity can be written 



E n (e-^-^ -i)-E n ^^'') h^o{{s.^)^^^^ - . 

Gee. {ijjeG GeT„ {i,j}eG 

(6.19) 

The full definition of the measure \q can be found in [3]; here we only mention 
its relevant properties. \q depends on the tree G and it is a probability measure 
supported on a subset L C {s^ G [0, 1] : 1 ^ i < 3 ^ n\. If satisfies the 
stability condition, then 

E SyReuy ^ -^6,, (6.20) 

,i<j iGl 

for all {sij} G L and all / C {1, . . . , n}. 

If Re Uij oo for any we immediately get the claim by using the stability 
condition above, since / dXcdsij}) — 1. The extension to possibly infinite numbers 
can be obtained using a trick due to Procacci [22]. Let G be a fixed tree and e > 0. 
We introduce 

{(l-£)m if G GnP, 

e if{*,j}GG\P, (6.21) 

otherwise. 
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Then 

1 ^ z<j ^ n 1 ^ i<j ^ n 1 ^ i<ij ^ n 

(6.22) 

If m is large enough (depending on e), {w-J™^} is stable and the first term of the 
right side is bounded below by — ^i- Then 



{i,i}eG 



(m.e) £ \ /• 

n ""^^ ' Klj ydAG(K-})e-^-<.-^-'r'. (6.23) 



A special case of the tree identity (|6.19p is 

{z,j}eG {i,j}6G (6.24) 

We get 



n (-^ir^) / dAG({%})e 



{iJ}eG 



{»j}eG\p 



We can insert this estimate into Eq. (|6.19p . Letting m ^ oo and then e ^ 0, we 
get Proposition [nH] (b). □ 

7. Proofs of the theorems 

In this section we prove the theorems of Section [2l We consider only the case 
where Assumption [2] holds true — the case with Assumption [2l is entirely the same, 
one only needs to replace all \C{x,y)\ with \u{x,y)\ and all e^''^'^ with e^^'^ . The 
proofs are based on the following tree estimate, which is a direct consequence of 
Proposition 16. II for almost all xi, . . . , x„ g X, 



n 

|^(a;i,...,a;„)| ^ -He^'^^') E niC(^-^^-)l- (7.1) 
Proof of Theorem \2.1\ We start by proving the bound (|2.9p . Let us introduce 

^^(^i) = Er^TM / diH(^2)... /diMi(x„)ne^''(^'' E n ic(^-^.)i, 



1=1 ggt„ {ij}eG 



K{x) = lim Kn{x). (7.2) 

(The term n = 1 is equal to e^'''-^^-' by definition.) We show by induction that 

KNix) < e"(")+2''(") (7.3) 
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for any TV. Then K{x) < for almost all x, and using jTl]) we get jM]). 

The case = 1 reduces to 1 ^ e"'^' and it is clear. The sum over trees with 
n vertices can be written as a sum over forests on {2, . . . , n}, and a sum over edges 
between 1 and each tree of the forest. Explicitly, 



Kn{xi) = X! fyt„ 1)1 51 / dl^K^a)... / d|^|(x„) 

"=1 ''m^l {Vi,...,V^} 

partition of {2,...,n} 

]-[e2K-.) |C(:ri,x,)| n IC(^-^.-)l)- (7.4) 

Here, T{V) denotes the set of trees with V as the set of vertices. If \Vk\ = 1 the 
sum over G £ TiVk) is one by definition. The term after the sum over partitions 
depends on the cardinalities of the V^'s, but not on the actual labeling. Also, each 
i &Vk gives the same contribution. We get 



n— 1 m ^ 1 ni,...,7im ^ 1 k—1 ^ ^ ' 
\'njri,—rL—\ 

/dlA^Kyi)... /d|A.|(y„J|C(xi,yi)ine2''(^') 11 IC(y.'%)l) (7-5) 

We obtain an upper bound by releasing the constraint ni + . . . rim ^ — 1 to 
fik ^ iV— 1,1 ^ fc ^ TO. We then get 

K^(x^) ^ e2^(-^)exp|^^^-^ ^ d^yi) . . . J d|/i|(y„) |C(a:i, yi)| 

f[,2Ky.) J2 n \ay.,y,)\] (7-6) 

i=i Ger„ {ij}eG ^ 

- e^''^-^) exp( / d\^l\{y,)\axuyl)\KN^,{y,) 



We have ifjv-i(yi) ^ e^'^^^+^^^^i^ by the induction hypothesis. Eg. (fTSl) follows 
from Assumption [2l 

The rest of the proof is standard combinatorics. The partition function can be 
expanded so as to recognize the exponential of connected graphs. Namely, we start 
with 

^=1+ E ^ /d^(^i)--- /'^^^''"^ ^ n C(x.,x,). (7.7) 

The graph G & Gn can be decomposed into k connected graphs whose sets of vertices 
form a partition of {1, . . . ,n}. Summing first over the number rrii of vertices for 
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each set of the partition, we get 



k\ ^-^ mil . . . m]^\ 

^ 1 k ^ 1 mi,...,m/j^l 
mi-\ {-ni/^—n 



= 1 



k 

EE^ E n{ /"d/x(a;i)... /'d/i(x„J(p(a;i,...,x™j|. 

n^lfe^l ' mi,...,mfc ^ 1 £=1 
JTiiH \-7nk—n 

(7.8) 

The triple sum is absolutely convergent thanks to the estimate p.9p that we have 
just established. One can then interchange the sums by the dominated convergence 
theorem. This removes the sum over n, and this completes the proof of Theorem 
O □ 



Next we prove Theorems 12.21 and 12.31 in reverse order, since we will use the 
convergence properties in the latter theorem to get the former. 

Proof of Theorem \2.3[ From the definition (|2.13p and the tree estimateEU we have 

\Z{x,,x,)\ ^ 7^^, j d|A*l(x3)... / d|Ai|(x„)ne2^(-') ^ [] KC^^-^^OI 

(7.9) 

The expression above involves a sum over trees of arbitrary size that connect 1 
and 2. Any such tree decomposes into a line of m + 1 edges that connect 1 and 2 
(m ^ 0), and m + 2 trees rooted in the vertices of the connecting line. Taking 
into account the combinatorial factors, we obtain 

\Zix,y)\ < \C{x,y)\K{x)K{y) 

+ J2 fd\fi\{x,)... fd\ti\{x„,)(fl\ax^,x,+i)\K{xi))K{y) 

with xq = X and Xm+i = y- The result follows from the bound (j7.3p for K . □ 

Proof of Theorem \2.Sl It is actually similar to the end of the proof of Theorem 
12.11 Z{x) can be expanded as a sum over graphs, that can be decomposed into a 
connected graph that contains 1, and other connected graphs. Taking into account 
the combinatorial factors, the contribution of connected graphs containing 1 yields 
Z{x), and the contribution of the others yields the expression (|2.8p for Z. Thus 
Z{x) = Z{x) Z. One step involved interchanging unbounded sums, which is justified 
because everything is absolutely convergent, thanks to (|2.9p and Theorem 12.31 

In the graph expansion for Z{x,y), the terms where 1 and 2 belong to the same 
connected graph yield Z(x, y) Z , and the terms where 1 and 2 belong to different 
connected graphs yield Z{x)Z{y)Z . The detailed argument is the same as above. 
We then obtained the desired expression. □ 
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